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ABSTRACT 

We model the expected X-ray polarisation induced by complex reprocessing in the 
active nucleus of the Seyfert-2 galaxy NGC 1068. Recent analysis of infrared interfer- 
ometry observations suggests that the ionised outflows ejected by the central engine 
are not aligned with the symmetry axis of the obscuring torus. This conclusion was 
obtained by extrapolating the apparent orientation of the narrow line region to the in- 
ner parts of the ionisation cones. We show that future measurements of the soft X-ray 
polarisation vector unambiguously determine the orientation of the ionisation cones. 
Furthermore, X-ray polarimetry across a broad photon energy range may indepen- 
dently verify the misalignment between the ionisation cones and the axis of the torus. 
To model the expected polarisation percentage and position angle, we apply the radia- 
tive transfer code stokes. Reprocessing of the primary X-ray radiation takes place in 
the accretion disc, the surrounding equatorial torus and the inclined, ionised outflows. 
We also examine additional equatorial scattering occurring in between the accretion 
disc and the inner surfaces of the torus. Radiative coupling between the different re- 
processing components is computed coherently. The resulting polarisation properties 
depend on the optical depth of the reprocessing regions and on the viewing angle 
of the observer. We show that even under unfavourable conditions the misalignment 
of the outflows with respect to the torus axis can be determined from a rotation of 
the polarisation position angle between softer and harder X-rays. We argue that the 
misalignment of the outflows with respect to the torus axis in NGC 1068 may be 
constrained by a future X-ray mission if equipped with a broad band polarimeter. 

Key words: radiative transfer - polarisation ~ galaxies: Seyfert - galaxies: individual: 
NGC 1068 - X-rays: galaxies - techniques: polarimetric 



1 INTRODUCTION 

The standard unified scheme for active galactic nuclei 
(AGN) attempts to explain the observational appearance 
of different AGN types mainly as an inclination effect 
l|Antonuccil [l993l : lUrrv fc Padovanilll995l ). The continuum 
emission being produced by the accretion disc of the central 
supermassive black hole (SMBH) irradiates the surrounding 
broad line region (BLR) and gives rise to broad optical/UV 
line emission. Obscuring equatorial dust at larger distances 
plays a crucial role in the unified model as at higher incli- 
nations it covers the BLR and thereby disentangles type- 
1 AGN, which spectroscopically show broad line emission, 
from type-2 objects, which do not. The accretion process 
onto the SMBH partly leads to re-ejection of matter and 
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supposedly causes ionised winds in the polar direction. The 
presence of these winds has brought strong support to the 
unified model as it allowed to indirectly detect hidden BLRs 
in obscure d type-2 objects by scatter ing of BLR light around 
the torus (|Antonucci fc Milled 1 1985l V Systematic studies in 
optical spectropolarimetry lead to a classification of Seyfert 
galaxies accor ding t o the direct ion of the opt ical polarisation 
angle (see e.g. lAntonucci 198j: [Smith et al.l i2002. 2004. and 
references therein) that can be explained in the framework 
of the unified scheme when including equatorial scattering 
inside the torus funnel. 

The simplest approach of the unified model assumes 
that the rotation axis of the accretion disc is aligned to 
the axis of the torus and of the outflows. The alignment 
of these components is justified by assuming a symmetric 
mass transfer from the inner torus towards the outer ac- 
cretion disc as well as by the symmetric coUimation effect 
that the torus should have on the polar outflows. But recent 
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work by iRaban et all (|2009l ) surprisingly suggests that the 
outflows of the well-studied Seyfert-2 galaxy NGC 1068 are 
inclined by 18 ° with respec t to the apparent axis of the torus 
(see Fig. 9 in lRaban et al.|[200S l. If this result is confirmed 
it has important consequences for our understanding of the 
accretion and ejection mechanisms close to active SMBHs. 

In this work, we explore to which extend the misalign- 
ment of the outflows with respect to the torus in NGC 1068 
can be probed by X-ray polarimetry observations. With the 
advent of a new generation of X-ray polarimeters based on 
photo-electric effects, X-ray imaging polarim etry becomes 
feasible across the 2-35 keV energy range (|Soffltta et al.l 
I2OIOD . The X-ray emission from AGN is much less sensi- 
tive to stellar contributions from the host galaxy than op- 
tical/UV radiation is, and therefore it offers a better probe 
of the reprocessing geometry in the innermost regions of 
NGC 1068. One could expect that the radiative coupling 
of centrally emitted X-rays to the torus and the outflows 
is somewhat similar to the observed behaviour in the opti- 
cal/UV as in both wavebands electron scattering plays an 
important role. But in the hard X-ray band electron scat- 
tering becomes wavelength-dependent, which should have an 
impact on the polarisation. 

In the following, we build up a sophisticated X-ray re- 
processing model of NGC 1068. We start out with a simple 
irradiated slab representing the accretion disc and an ele- 
vated primary source. Then, we add dust obscuration in a 
torus geometry testing for different optical depths. Next, we 
include inclined polar outflows and conduct calculations for 
different optical depths of these ionisation cones. In the soft 
X-ray band with photon energies E < 10 keV absorption 
processes in the torus are very important and at high in- 
clination angles the primary source and reflection from the 
accretion disc should only be seen indirectly by scattering in 
the outflows. The scattering in the ionisation cones should 
then imprint a polarisation position angle that is perpen- 
dicular to the axis of the outflow. At higher energies, when 
scattering off the torus becomes more prominent, the posi- 
tion angle should change and be determined by the torus 
geometry. This switch of the polarisation angle depends on 
the relative Stokes fluxes coming from the ionisation cones 
and from the torus. Therefore, it must vary with the opti- 
cal depths in both components as well as with the inclina- 
tion angle. Finally, we also examine the impact of additional 
equatorial scattering in the farther away components of the 
accretion flow that are located between the inner surfaces of 
the torus and the outer accretion disk. Such an additional 
scattering component produces polarised reprocessing by it- 
self but also changes the irradiation geometry of the inner 
surfaces of the torus. 

We test the resulting spectra and polarisation proper- 
ties for different optical depths of the obscuring dust and of 
the scattering electrons in the ionisation cones. Aside from 
the optical depths, our model comprises other interesting 
parameters such as the half-opening angle of the torus or 
of the outflows. Their impact is going to be examined else- 
where. This study is particularly designed to predict the 
ex pected X-ray polari sation signal if the geometry inferred 



scribed in Sect. [S] The results for different combinations of 
the individual reprocessing regions and then for the case of 
NGC 1068 follow in Sect. [J] We close with a discussion in 
Sect. [5]where we also shed light on the prospects for future 
X-ray polarimetry observations of NGC 1068. 



2 X-RAY MODELLING WITH STOKES 

We model the radiative transfer and the polarisation signa- 
ture of the complex X-ray reprocessing environment in AGN 
by applying the code STOKES ( Goosmann & GaskcU 2007). 
It allows to define rather general emission and reprocessing 
regions and coherently computes the 3D radiative coupling 
between them. The model space can then be evaluated from 
any viewing direction as a function of the photon energy, 
E, in terms of the X-ray spectral flux, F/F,, the polarisa- 
tion degree, P, and the polarisation position angle, ip- Note 
that the flux is normalised to the total flux of the primary 
source. Ft, emitted into the same viewing direction. The 
STOKES code was originally developed for the optical/UV 
spectral r ange and has rece ntly been extended into the X- 
ray band l|Goosmannl I2OO9I ) . We generalised the algorithm 
for electron scattering to include the inelastic Compton ef- 
fect and we included routines to treat photo-ionisation and 
recombination effects. Photo-absorption above the iron K- 
shell and the subsequent emission of iron Ka or K/3 line 
photons is included and weighted against the probability of 
Auger effects. Throughout this paper, we consider X-ray re- 
processing by neutral matter consisting of hydrogen, helium, 
carbon, nitrogen, oxygen, neon, magnesium, silicate, sulfur 
and iron with cosmic abundances. 

The X-ray part of the code was tested against pre- 
vious results of spectral modelling, in particular the re- 
cent X-ray reproce s sing m odel for AGN tori presented by 
Murphv fc Yaaoog (120091 ) and disc reflection models by 
George fc FabianI (|l991f ) . Polarisation results were compared 
against the results given by Ghisellini, Haardt & Matt 
(1994). 



bv lRaban et al.l l|2009l ) is correct. We therefore fix all dimen- 
sions and angles according to their findings. 

The radiative transfer code applied here is presented 
in Sect. [2] and the details of the model geometry are de- 



3 X-RAY REPROCESSING IN NGC 1068 

The modelling geometry we consider is based on several el- 
ements. We approximate the reprocessing of an AGN accre- 
tion disc jGeorge fc Fab"ianlll99ll : iMatt. Perola fc Piro|[l99ll : 
lMattlll99a r by defining a flat, cylindrical, optically thick 
and neutral disc that is being irradiated from both sides by 
point-like X-ray sources isotropically emitting a power-law 
spectrum with a spectral flux Ft{E) oc E~°' and a = 1. 
Note that this spectral index a corresponds to a photon in- 
dex r = 2. This choice of the spectral index is motivated 
by the analy sis of an XMM-N ewton spectrum of NGC 1068 
presented in lMatt etlii] (|2004l '). The two primary sources in 
our model are located on opposite sides of the disc symme- 
try axis such that the disc subtends a half-opening angle 
of ~ 75° with the sources. The reprocessing medium thus 
catches more than 95 per cent of the primary radiation emit- 
ted towards the disc. 

To this base element of our modelling scenario we then 
add up to three more reprocessing regions: an equatorial, ob- 
scuring torus with an elliptically-shaped cross-section (1), 
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Table 1. Characteristics and dimensions of the reprocessing elements used in our modelling. 
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reprocessing region 


description 
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and half-height of 3.25 X 10~^ pc. 




irradiated by two isotropic sources at 




iO.OOl pc from the disc mid-plane. 




filled with optically-thick, neutral matter 


equatorial, obscuring torus 


torus with elliptical cross-section and half-opening angle 9tor = 60° , 




inner-radius of 0.1 pc and outer radius of 0.5 pc, 




filled with neutral matter of an equatorial, hydrogen-equivalent density A^h tor 


tilted, conical outflows 


double-cone with the half-opening angle Scone = 40° , 




titled by the angle ^tin = 18° from the symmetry axis. 




inner-radius of 0.3 pc and outer radius of 1.8 pc. 




filled with electrons of a radial Thomson optical depth Tconc 


equatorial, scattering disc 


horizontal, flared disc with half-opening angle (^wodgc = 20° , 




with ^wcdge being measured from the equatorial plane. 




inner-radius of 0.03 pc and outer radius of 0.05 pc. 




filled with electrons of a radial Thomson optical depth r^cdgc = 1 



tilted, conical outflows in the polar directions (2), and a 
ring-like scattering region in the equatorial plane, located 
between the accretion disk and the inner surfaces of the 
torus (3). The equatorial torus is filled with neutral matter 
and centred on the much smaller reprocessing disc. The di- 
mensions of all emission and scattering geometries are sum- 
marised in Table [T] and a sketch is provided in Fig. [1] The 
symmetry axes of the torus, the equatorial wedge and the 
disc are identical and the tor us funnel has a half- opening an- 
gle 6'tor = 60°, as taken from lRaban et al] |200^. Thus, the 
direct view of the central components is obscured at viewing 
directions, i, with i > 6'tor. Both angles are measured with 
respect to the symmetry axis. 

We first only consider reprocessing off the irradiated 
accretion disc alone and then add the torus. The following 
scenarios include polar outflows that we realise by a tilted 
double-cone. The centre of the cone is identical with the 
ones of the reprocessing disc and the torus. The double- 
cone half-opening angle is measured by Scone = 40° and 
its tilting angle w ith respect to the torus symmetry axis 
by Stilt = 18° (see iRaban et al.l l2009l'). We assume that the 
matter inside the double-cone is ionised and dominated by 
electron scattering. The Thomson optical depth, Tconc, is 
measured between the inner and outer radii of the cone. The 
misalignment between the torus and the ionisation cones is 
inferred from the m id-IR interferometry measurements of 
iRaban et all (|2009|) in combination with the maser obser- 
vations of Gallimore et al.l (|l996l) as well as with the HST- 
spectroscopy modelling conducted bv lDas et all (|2006l ). For 
further details on the geometrical consi derations we refer t o 
the very instructive discussion given in iRaban et all (|200i). 

There is another, potentially important scattering re- 
gion that is required to explain the dichotomy of the opti- 
cal/UV polarisati on ang l e. Th e equ atorial scatte r ing ri ng as 
described by Ant onuccif (Il984l ') and lSmith et all (|2004l 'l effi- 
ciently polarises with a projected polarisation vector being 
parallel to the symmetry axis. The ring should also con- 
tribute in the X-ray range and we model it as an equatorial 
scattering wedge between 0.03 and 0.05 pc. The symmetry 
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Figure 1. Sketch of the emission and scattering geometry in this 
modelling. The insets show a consecutive zoom on the central 
elements representing a neutral slab that is surrounded by an 
equatorial scattering wedge and irradiated on both sides by point- 
like, primary X-ray sources. 



axis of the wedge is aligned with the torus axis and its half- 
opening angle with respect to the equatorial plane is 20°. 
The equatorial scattering region is assumed to be at least 
partly ionised and we therefore assume pure electron scat- 
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Figure 2. Modelling results for the reprocessed spectrum of an ir- 
radiated, neutral accretion disc. Top: spectral flux, F, normalised 
to the centrally emitted flux, as a function of the disc incli- 
nation i measured by /z = cosi. Bottom: polarisation degree, P. 
The polarisation vector, i/i, is aligned with the projected disc axis 
at all inclinations (i/i = 90°). 



tering with a Thomson optical depth of unity between the 
inner and outer radius of the wedge. 

When showing modelling results we take into account 
existing (anti-)symmetries between the two hemispheres 
above and below the equatorial plane. We therefore only 
discuss the results for viewing directions 0° < i < 90° (or 
1 > cosi > 0) with i being measured with respect to the 
symmetry axis of the disk and the torus. 
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Figure 3. Modelling results for a system of an irradiated accre- 
tion disc and an equatorial torus (see Table [TJ with N^i tor = 
lO'^'^cm"^. From top to bottom the flgure shows the spectral flux 
in terms of F/Ft, the polarisation percentage P and the polari- 
sation position angle -0 as a function of the observer's inclination. 



4 X-RAY POLARISATION RESULTS 

4.1 The reprocessing spectrum of the irradiated 
accretion disc 

The spectral modelling of the accretion disc irradiated by an 
elevated X-ray source confirms previous results, as is shown 
in Fig. [21 The reprocessed spectra reveal iron Kq and K/? 
fluorescence lines, the associated iron K absorption edge and 
the broad Comptonised hump centred around 30 keV. Ab- 
sorption efi^ects mainly due to hydrogen and helium produce 
an overall positive slope of the normalised reprocessed spec- 
trum. Note that in the figure we have omitted the contribu- 
tion of the directly visible, unpolarised primary radiation to 
have spectral features in the flux and polarisation spectrum 
come out more clearly. 

In this modelling case, the polarisation vector, ^j), is al- 
ways aligned with the projected symmetry axis. At face-on 
view, when i is low, P does not exceed 1 per cent because 
the scattering geometry is almost symmetric. The polarisa- 
tion degree then rises with increasing i and the scattering 
medium appears less symmetric with respect to the line- 
of-sight. The polarisation degree drops sharply across the 
iron K lines, which is due to dilution by the unpolarised flu- 
orescent emission. Across the Compton hump the relation 



between P and i differs from the soft X-ray band, which is 
due to the Compton scattering phase function that favours 
forward over backward scattering at higher photon energies. 

It is instructive to co mpare these results to the non- 
relativistic calculations bv iMatj (|l993l ). Qualitatively, we 
find the same behaviour of the polarisation degree and an- 
gle but the absolute values of P obtained h ere ar e lowe r 
than for comparable cases shown in fig. 4 of iMattI (|l993t ). 
This can be explained when considering that the net po- 
larisation degree of the reprocessed radiation results from 
integrating the reprocessed Stokes fiux over the whole disc. 
After one scattering event, the polarisation of the outgoing 
radiation is directed perpendicularly to the scattering plane, 
and thus the central parts of the disc produce polarisation 
angles around tj) = Q° (perpendicular to the symmetry axis) 
while the outer regions of the disc rather give rise to t/) = 90° 
(parallel to the projected symmetry axis). The net Stokes 
fiux is dominated by the outer parts of the disc but still in- 
fiuenced also by the perpen dicular polar isation state coming 
from the disc centre. Since iMattI ljl993h included a central 
hole in the disc, the impact of the disc centre on the po- 
larisation was smaller and the net polarisation percentages 
found were larger. 

Detailed calculations of the reprocessing by an irradi- 
ated accretion disk may also consider the presence of the 
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marginally stable orbit and relativistic ray-tracing effects 
around th e black hole. Such resu l ts hav e rec ently been pre- 
sented by ISchnittman fc KrolikI ll201Cf) and Dovciak et alj 
(|201ll ). The treatment bv lSchnittman fc KrolikI l|2010h even 
involves secondary reprocessing, as the reflected component 
can be directed back to the accretion disc by relativistic 
light-bending. Note that in the work presented here we as- 
sume that the source is located high enough above the disc 
so that the lack of reprocessing due to the central hole as 
well as relativistic effects can be neglected. 



4.2 Reprocessing by the torus 

The resulting normalised flux and polarisation spectra for 
the combined reprocessing by an irradiated accretion disc 
and an equatorial obscuring torus are shown in Figs. 3-4 
for a torus hydrogen column density A^n.tor = lO^^cm"^, 
A'n.tor = lO'^^cm"^ and A^n.tor = lO^^cm"^, respectively. 
The flux spectra reproduce the features of neutral repro- 
cessing but differ signiflcantly above and below the torus 
horizon: for cosi > 0.5 the irradiating primary source dom- 
inates the resulting spectrum. The underlying power-law 
with a = 1 appears as a horizontal curve in the F/ Ft repre- 
sentation and the additional reprocessing features are rela- 
tively weak. For cos i < 0.5, the primary source is no longer 
directly visible and soft X-ray absorption becomes more pro- 
nounced as A^'n.tor increases. The amount of absorption in- 



creases with i as the radiation must undergo more scattering 
events to escape from the torus funnel toward high inclina- 
tions. 

The behaviour of the polarisation degree and an- 
gle varies strongly w ith A'^e.tor. As was pointed out by 
iGhisellini et al.l (| 19941 ) for a different torus geometry, the net 
Stokes flux is determined by both reflected and transmitted 
components that differ in their direction of the polarisation 
angle. For the case of A'n.tor = lO^^cm"^ absorption effects 
are rather small and mostly apparent in the soft X-ray spec- 
trum at the highest inclinations. Otherwise, the polarisation 
is determined by transmitted radiation penetrating through 
the torus walls. The intrinsic polarisation with ij) = 90° com- 
ing from the accretion disc is thereby reinforced because of 
scattering in the rather flat geometry of the torus, which has 
a considerably large half-opening angle of 60° . 

When absorption effects become important, then the re- 
processed radiation is more dominated by scattering off the 
torus' inner walls. Since the absorption cross-sections are a 
strong function of photon energy, the polarisation position 
angle can switch between the two orthogonal states -0 = 0° 
and = 90°. Beyond ^ 12 keV, where the integrated ab- 
sorption effects become less important than the efficiency of 
electron scattering, the hard X-ray polarisation angle tends 
always to be parallel to the projected torus axis. When going 
through the switch in polarisation angle, P is near to zero. 
The polarisation degree also diminishes across the iron K 
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lines due to strong dilution by the intrinsically unpolarised 
line flux. 

The polarisation degree and position angle depend 
strongly on the inclination. For a pole-on view, the system 
of the irradiated disc and torus is close to being axisymmet- 
ric and the resulting polarisation of the reprocessed radia- 
tion is low. Additionally, the direct unpolarised flux from 
the primary source further lowers the observed polarisation. 
When the torus is sufficiently optically thick, as in the case 
of Nn^tor = lO^'^cm"^, the net polarisation state is mainly 
determined by scattering off the inner surfaces of the torus 
funnel. For a given line-of-sight towards the torus, the rela- 
tive sizes of its visible inner surfaces oriented perpendicular 
and parallel to the line-of-sight determine the resulting po- 
larisation angle. Scattering off surfaces which are parallel to 
the line-of-sight induces aligned polarisation with tp = 90°, 
the perpendicular surfaces, on the other hand, produce per- 
pendicular polarisation with tp — 0° . When increasing the 
viewing angle i, the reflected radiation coming from per- 
pendicular surfaces is mainly due to forward and backward 
scattering, which leads to a weaker Stokes flux for the per- 
pendicular polarisation component. At the same time, the 
relative exposure of parallel and perpendicular surfaces with 
respect to the line-of-sight also changes with i. The interplay 
between both effects can be complicated as is witnessed by 
the partly complex spectral shapes of t/i shown in Figs. 3-4. 
The explanations given here are similar but not equivalent 
to the modelling results o btain ed in the optical/UV band 
and d iscussed by iKartid (|l995h and |g oosmann fc Gaskelll 
l|2007l ). The differences to the discussion of the optical/UV 
range are a result of the particular energy dependence of 
the absorption and scattering cross-sections as well as of the 
scattering phase functions in the X-ray band. 

4.3 Reprocessing by the disc, a thick equatorial 
torus and polar cones 

To include the effects of polar scattering, we add double- 
conical electron scattering regions with low optical depth 
to the previous model of an irradiated accretion disc and 
equatorial torus. We evaluate the most opaque of our torus 
models. This choice is motivated by the X-ray spectroscopy 
data of NGC 1068 confirming that the hydrogen column 
density deduced from the 0.2-10 keV band must be larger 
than lO^^cm"^ (Matt et al. 1997; Bianchi, Matt & Iwasawa 
2001). Depending on the observer's inclination, the hydrogen 
column density integrated along the line-of-sight can be by 
a factor of a few lower than the column density A^H.tor that 
is defined along the equatorial plane. By setting A'^e.tor = 
lO^^cm"^ we make sure that the torus remains sufficiently 
optically thick at all inclinations above 6toi = 60°. 

We take into account that the double-cone is inclined 
by 18° with respect to the torus symmetry axis. We assume 
a Cartesian coordinate system with the torus symmetry axis 
being the 2-axis and the xy-plane being the torus equatorial 
plane. Then, the double-cone symmetry axis is tilted inside 
the yz-p\a,ne leaning toward the positive y-axis. The angle 
Scone is measured between the double-cone axis and the z- 
axis. Two different Thomson optical depths, Tconc = 0.03 
and Tcone = 0.3 are considered. 

The results are plotted in Fig. (5] For each inclination, we 
consider the azimuthal viewing angle, 0, along the positive 



i-axis integrating the flux over 20° in azimuth. Such a view 
corresponds to the illustration given in Fig. [T] Compared 
to the models without polar scattering, the soft X-ray spec- 
trum seen at i > 6'tor is much less absorbed. Even at low op- 
tical depth, the double-cone scatters an important fraction 
of the primary emission around the torus towards high incli- 
nations. This is the 'periscope effect' of polar scattering that 
has been evalu ated and observa tionally expl oited in the opt i- 
cal/UV range (|Antonucci fc Miller 1985 ; Smith et al.ll2002l '). 
Since the periscope effect is due to polar electron scattering, 
it is not surprising that it also exists in the X-ray range. 
Even at low Tcone the continuum polarisation vector in the 
soft X-ray band is dominated by polar scattering inducing 
a polarisation angle oi tp = 18° that corresponds to the in- 
clination of the outflows. The soft X-ray polarisation rises 
with i reaching a maximum of about 60 per cent at edge-on 
view. 

At higher photon energies, the polarisation induced by 
the torus becomes more important because its albedo rises. 
The polarisation spectrum is again highly influenced by di- 
lution across the iron K lines, but past the iron K absorption 
edge the Compton scattering off the torus becomes dominant 
for intermediate inclinations. The exact outcome of the po- 
larisation depends on the relative Stokes fluxes coming from 
the polar outflows, on the one hand side, and from the ac- 
cretion disc and torus, on the other hand side. When viewed 
strictly edge-on the absorption effects in the torus are too 
strong and the polarisation angle is always determined by 
the polar outflows. But a rotation of the polarisation an- 
gle between the soft and the hard X-ray band is observed 
and could be used to measure the misalignment between the 
symmetry axes of the torus and the outflows. In the flrst and 
second column of Tabled we quantify the difference in 
the polarisation angles between 20 keV and 2 keV for differ- 
ent obscured viewing angles. The values of Atp vary strongly 
with the inclination but even at edge-on view and for a high 
optical depth of the ionisation cone a small rotation of the 
polarisation angle is noted. 

For contrast, we show in Fig. [5] the results for the X- 
ray polarisation of a model including polar outflows with 
Tcone ~ 0.03 but being aligned to the 2-axis. Without in- 
clining the double-cone, tp can only adopt the two possi- 
ble values that correspond to the perpendicular polarisation 
caused by the polar scattering in the soft X-ray band and the 
aligned polarisation vector due to the disc/torus at higher 
photon energies. This bimodal distribution of the polarisa- 
tion vector is related to the symmetry of the model with 
respect to the z-axis. The comparison shows that the grad- 
ual change of tp with photon energy, as it is apparent in the 
modelling case of Fig. \5\ (left), is characteristic for a mis- 
alignment of the outflows with respect to the torus. 

4.4 The effect of equatorial scattering 

The third and fourth column of Table [2] show some exem- 
plary results for a scenario that includes equatorial scatter- 
ing. The corresponding models are based on the cases of 
Sect. 14.31 but additionally include an equatorial wedge as 
described in Table [1] The additional scattering region has 
an impact at all viewing directions. For polar outflows with 
low optical depth Tconc, the rotation Atp of the polarisation 
angle can either rise or fall with respect to the models of 
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Figure 5. Modelling results for a system of an irradiated accretion disc, an equatorial torus with Affj.tor = lO'^^cm"'^ and inclined polar 
electron scattering cones with Tcono = 0.03 (left) and Tcono = 0.3 (right). See Table[T]for model details. The denomination of the panels 
is as in Fig. \3\ 



Table 2. Rotation Aip = V(20 keV) — ip{2 keV) of the polarisation angle and range of the polarisation percentage P between 20 keV 
and 2 kcV for the two values of Tcono examined in Sect. (column 1 and 2) and Sect. H7H (column 3 and 4). 



viewing angle i 



Tcono — 0.03 

Alp, AP 



: 0.3 



: 0.03 (cquat. wedge) 



: 0.3 (equat. wei 



Atp, AP 



Ail, AP 



Aip, AP 



: 0.45, 



: 0.35, 



: 0.25, 



: 0.15, 



: 0.05, 



; 63° 



; 70° 



; 76° 



; 81° 



■ 87° 



84.9° - 20.0° : 
4.8% - 21.3% 
74.6° - 18.3° : 
2.8% - 41.3% 

40.6° - 17.8° : 
4.0% - 50.3% 



: 64.9° 



: 56.3° 



: 22.8° 



23.4° - 
16.6% 



17.5° = 5.9° 
- 56.9% 



18.7° - 17.5° = 1.2° 
52.7% - 60.2% 



32.9° - 19.3° = 
6.8% - 36.7% 
23.4° - 18.4° = 
13.6% - 45.6% 

20.3° - 17.9° = 
25.8% - 51.0% 

18.6° - 17.6° = 
43.5% - 55.1% 

17.9° - 17.3° = 
55.3% - 57.7% 



13.6° 85.9° - 78.4° = 7.5° 
7.7% - 6.9% 

5.0° 75.8° - 20.1° = 55.7° 
4.5% - 41.3% 

2.4° 47.8° - 19.3° = 28.5° 
4.4% - 48.8% 

1.0° 24.8° - 18.8° = 6.0° 
17.0% - 55.2% 

0.6° 18.8° - 18.5° = 0.3° 
53.1% - 58.8% 



39.9° - 27.3° = 12.6° 
7.2% - 15.9% 

24.7° - 20.0° = 4.7° 
15.8% - 44.6% 

20.9° - 19.4° = 1.5° 
26.4% - 49.6% 

19.2° - 18.9° = 0.3° 
43.8% - 53.7% 

18.4° - 18.6° = -0.2° 
55.3% - 56.4% 



Sect. 14.31 for a high Tcone, the value of Atj) always diminishes. 
At edge-on inclination, the polarised flux escaping from the 
torus funnel is too small to efficiently counterbalance the ef- 
fect of polar scattering. Therefore, the resulting A^/) is very 
low. 

The results illustrate that the wavelength-dependence 
of the scattering efficiency and the radiative coupling be- 
tween the different reprocessing regions are important to un- 
derstand the net polarisation. The equatorial electron scat- 



tering competes with the polar electron scattering in the 
outflows, independently of spectral energy. Therefore, the 
polarisation position angles are systematically higher when 
equatorial scattering is included. The rise in il) is more ap- 
parent at intermediate inclinations but still exists towards 
an edge-on view. 

The additional scattering region also has an impact on 
the rotation A'i/) between the soft and hard X-ray polari- 
sation angle. The equatorial wedge covers a fraction of the 
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Figure 6. Polarisation percentage and position angle as a func- 
tion of photon energy for a model of an irradiated accretion disc, 
an equatorial torus with Afji.tor = lO^^cm"^ and polar scatter- 
ing cones with Tcono = 0.03 seen at the obscured viewing angles 
fi = 0.45 and ^ = 0.25. In this case, the polar outflows are aligned 
to the z-axis and the results can be contrasted to the case of 
Fig. [5] (left). 



primary radiation going towards the inner surfaces of the 
torus and it changes the irradiation geometry for both the 
torus and the polar cones. These effects diminish the result- 
ing Ai/>, at least when the optical depth in the polar outflows 
is high. We have investigated different half-opening angles 
of the pure electron scattering wedge as well as a reduction 
of its optical depth. The obtained values of Ai/; change for 
such cases, but the qualitative behaviour stays the same. 



5 DISCUSSION 

5.1 Directly measuring the orientation of the 
inner outflows 

We have conducted accurate modelling of the expected X- 
ray polarisation induced by complex reprocessing in the 
active nucleus of NGC 1068. This work is motivated by 
the apparent misalignment of the ioni sation cones with re - 
spect to the torus axis as discussed bv lRaban et~al] (|2009D . 
The orientation of the ionisation cones was inferred from 
slit spectroscopy of the narrow line region (NLR) using the 
HST STIS in combination with kinetic modelling (iDas et al.l 
|2006| ). But|Raban ct al. (2009) also point out that this ori- 
entation is still a matter of debate. From the observation 
of infrared emission lines, Poncelet, Sol & Perrin (2008) de- 
duce a different position angle of the NLR that is more in 
agreement with the HST ima ging results of [OIII] emission 
reported by lEvans et al.l 11991). But then the UV-imaging 
probably suffers from significant foreground absorption. In 
addition to that, we want to point out that all these mea- 
surements are taken on large spatial scales reaching out to 
100 parsec from the central engine. It is not straightforward 



to derive from these observations what the actual geome- 
try of the ionisation cones across a few parsec is. One could 
imagine that the outflows become deflected or twisted on 
intermediate scal es and thus their o rientation would vary 
with distance. As iRaban et all (|2009h conclude, further in- 
vestigation is needed to really constrain the geometry of the 
innermost outflows. 

In this context. X-ray polarimetry is going to give un- 
ambiguous constraints. The scattered X-ray emission must 
come from very close to the central engine, as is indicated 
from X-ray spectroscopy results reve aling the presence of 
high-ionisation e mission lines of iron (iMarshah et al.lll993l : 
iMatt et ahlbooj ). Furthermore, any foreground absorption 
effects seen in the UV should much less interfere in the X- 
ray range so that X-ray polarimetry enables a direct view 
on the innermost parts of the outflow. Our modelling results 
have shown that the polarisation position angle in the soft 
X-ray range is clearly dominated by the polar scattering, 
which enables a precise measurement of the orientation of 
the scattering material. 



5.2 Determining the misalignment with the torus 

Another aim of this investigation is to verify if the rotation 
Ai/j of the X-ray polarisation angle between soft and hard 
X-ray energies can be used to measure the misalignment be- 
tween the polar scattering regions and the obscuring torus. 
We have found that coherent modelling of the different scat- 
tering components leads to values of A-i/) that at heavily ob- 
scured viewing directions allows one to distinguish the torus 
from the ionisation cones. The relation between Ai/) and the 
actual misalignment also depends on i. The highest values 
of A-(/i, which are most favourable for future X-ray polarime- 
try observations, are found at moderate viewing angles. At 
strictly edge-on viewing directions, Ai/i becomes very small, 
especially for higher optical depths of the ionisation cones. 
But, as we discuss in the following, our modelling approach 
is rather conservative, and in reality we should expect more 
favourable conditions for the measurement of /S.ip. The re- 
sults we present in this work give lower limits on At/i for 
NGC 1068 assuming the least favourable conditions for a ro- 
tation of the polarisation angle between low and high photon 
energies. 



5.2.1 Optical depth and geometry of the outflow 

The net Stokes flux and therefore the resulting polarisa- 
tion angle is determined by the competition between the 
disc/torus system and the polar scattering cones. Especially 
at higher inclination the polar scattering has a strong im- 
pact as the scattering-induced polarisation rises towards or- 
thogonal scattering angles. The efficiency of polar scatter- 
ing depends on Tcone and is low for optically thin scattering 
cones. It rises with increasing Tcone and then goes through 
a maximum until multiple scattering effects become impor- 
tant and diminish the resulting polarisation. Note also, that 
the Stokes flux from the double-cone results from integrating 
photons over its entire opening angle. This leads to combin- 
ing scattered photons having a large range of polarisation 
orientations and also diminishes the net polarised flux. Over- 
all, the case of Tcone = 0.3 that we investigate produces a 
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strong Stokes flux but at most inclinations it remains com- 
parable to the polarised flux coming from the disc/torus 
system and therefore the rotation Atp is significant. Having 
in mind the unified scheme of AGN, it is important that 
Tconc remains low enough as the ionised winds are seen in 
transmission in Seyfert-1 galaxies, where they produce warm 
absorption in the soft X-ray band. Observations of the warm 
absorber rule out that the medium has a hi gher column den- 
sity t han 10'^'^ cm^'^ for most objects (see [Turner fc Millerl 
I2OO9I . and references therein). Only a few exceptions exceed 
this threshold so that our choice of Tcone = 0.3, which corre- 
sponds to a column density around 4.5 x lO^^cm"^, is a good 
estimate for the upper limit. Note also that for NGC 1068 
the column density of the polar scattering medium has 
been constrained from the observed emission measure dis- 
tribution to 4 x lO^^cm" ^ < logiVH.Cono < 7 x 10^^ cm"^ 
iKinkhabwala et al.|[20o3 ). Furthermore, we assume in our 
modelling that the matter is uniformly distributed inside 
the double-cone, whereas the approach bv lDas et all (|2006l l 
is based on outflows with a shell-str ucture, such as predicte d 
for magnetically driven outflows ijKonigl fc Karti3 Il994h . 
For all these reasons, the Stokes flux coming from the ioni- 
sation cones of NGC 1068 should be lower than we assume 
in our modelling, which enables a larger A'i/). 

The misalignment of the ionisation cones wit h respect 
to the torus axis as derived bv lRaban et~all (|2009l ) refers to 
the projection of the double-cone onto the plane of the sky. 
But if the outflows are rotated in azimuth with respect to 
the j/z-plane, it is possible that the real misalignment with 
respect to the torus axis is larger than 18°. Note that a 
slight inclination of the double-cone towa rds the observer i s 
suggested by the kinematic modelling of iDas et al.l (|2006h . 
Then, the average scattering angle for the polar electron 
scattering is no longer close to 90° but will shift towards 
forward and backward scattering. This reduces the polarised 
flux coming from the outflows. At the same time, a stronger 
inclination of the cones also requires a larger half-opening 
of the torus, otherwise the outflows would be blocked. The 
scattering distribution of the torus is thus geometrically flat- 
ter in this case, which should increase its Stokes flux. Both 
effects work in favour of higher values for Atp than found 
from the models presented here. 

5.2.2 dumpiness of the torus material 

We are also pessimistic in our assumptions about the torus 
scattering properties. In our modelling we assume a max- 
imum optical depth by setting Nn,toT ~ lO'^^cm"'^. This 
agrees with previous spectroscopic studies showing that the 
obscuring material along the line-of-sight toward NGC 1068 
is entirely opaque. But for most viewing angles we are 
still more than one order above the lower threshold of 
the observed hydroge n column densi t y, wh ich was given by 
iMatt et al.l (|T997h and lBianchi et all (j200ll ) based on ASCA 
and BeppoSAX data. It is thus likely that the torus is less 
optically thick than we assumed in our modelling. A lower 
A'n.tor increases the polarised flux of the torus as we show 
in Fig.[7l Between the cases of lower and higher A^'n.tor there 
is an increase in Stokes flux at 20 keV by a factor of 2.5 at 
i « 76° and of more than a factor of 10 at i « 87°. 

The exact geometry of the obscuring dust region re- 
mains a matter of debate, although infrared interferom- 
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Figure 7. Modelled polarised flux spectra for a system of an 
irradiated accretion disc and an equatorial torus with TVjj tor = 
lO^^cm"^ (solid line) and A'jj.tor = lO^^cm"-^ (dashed line) for 
two different obscured viewing angles. 

etry observations start putting robust c onstraints o n the 
geometrical size of t he d usty region Jjaffe et al ] l2004l : 
iRamos Almeida et al.l l2009l : iHonig et all I2OI0I '). The latest 
generation of infrared radiative transfer models allows for a 
satisfying descrip tion of the data by as su ming that the dust 
region is clumpy (iNenkova et al.ll2008al lbl: ISchartmann et ahl 
l2008l : lHonig fc Kishimotdl2010l ). Such a clumpy composition 
of the torus could also address criticism that the unified 
model could not properly explain how the observed spatial 
extension of the dust region can be maintained against self- 
gravity. The clumpy structure would argue for a dynamical 
model in which the torus is not in hydrostatic equilibrium 
(jKrolik fc Begelmanlll988 l). In the results presented here we 
do not include the effect of dumpiness on the resulting polar- 
isation. But if the torus medium is clumpy and at the same 
time dynamical, there is a better chance for fluctuations in 
the column density along a given line-of-sight. Observing 
NGC 1068 with a lower opacity A^H.tor is thus more likely if 
the medium is clumpy instead of uniform. As shown above, 
a lower opacity enhances the Stokes flux from the torus and 
thus the rotation Atp. 

5.2.3 The primary source and equatorial scattering 

We assume that the polarisation coming from the accre- 
tion disc is based on reprocessing across the whole disk 
surface. The extended disk has a large radius compared to 
the height of the X-ray source above. Then, the outer disc 
regions produce polarisation with tp = 90°. This aligned 
polarisation is being diminished by the Stokes flux coming 
from the disk centre and having -i/i = 0°. As discussed in 
Sect. 14.11 such a setup assumes that the primary source is 
located far above the black hole so that the central hole of 
the disk does not have any significant impact. When relax- 
ing this assumption and considering sources that are much 
closer to the black hole, the reprocessing geometry changes 
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(see e.g. Martocchia & MattI 


19961: Miniutti & Fabianll2004 


ISchnittman & KrolikI 2010l: 


Dovciak et al.l 2011) and the 



missing part of the disc inside the marginally stable orbit 
becomes significant for the Stokes flux. The net polarisa- 
tion with i/) = 90° of the disc is thus stronger for such a 
case, which works in favour of a larger value for At/j. As- 
suming a different irradiation geometry, for instance a cen- 
tral but more extended hot corona, woul d not change this 
result much as is implied by the results of iMattI l| 19931 ). The 
main reason for this is that the disc remains very large with 
respect to the size of the corona. 

We have shown that generally equatorial scattering 
leads to a lower At/i, in particular at high inclinations of 
the observer. But note that we assume the equatorial scat- 
tering being entirely due to electrons. This is an extreme 
example as beyond the BLR also moderately ionised and 
neutral medium should exist. Less ionised matter has a dif- 
ferent scattering efficiency and should absorb more strongly 
in the soft X-ray band and thereby maintain higher values 
of A^. The two exemplary cases shown in Sect. [4^ rather 
define lower limits of the expected /S.xp when equatorial scat- 
tering is important. 



5.3 Encouraging prospects for future X-ray 
polarimetry of NGC 1068 

The observational window of X-ray polarimetry will soon 
be (re-)opened by the NASA Gravity and Extreme Mag- 
netism Small Explo rer ( GEMS) that i s currently prepared 
for launch in 2014 (|Swank et all l2009l : [jahoda 2010). The 
satellite will be equipped with a soft X-ray polarimeter 
reaching up to 10 keV. Confirmation of the polarisation po- 
sition angle induced by scattering in the misaligned outflows 
of NGC 1068 will thus principally be possible with GEMS. 
For the practical measurement, the collecting area of the 
mirrors might still be a limiting factor. On the other hand 
side, GEMS is entirely devoted to X-ray polarimetry and 
allows for long exposure times. 

The technology is also ready for mid-size observatories 
that include imaging polarimetry at high angular resolution 
and over the whole range of 2 keV to 35 keV. Such X-ray po- 
larimeters could efficiently observe near-by AGN. We use the 
response matrices of the polarimeters desi gned for the for- 
merly proposed New Har d X-ray Mission (jTagliaferri et all 
l20ld : ISoffitta et aLllioiol ). It turns out that a 300 ks ob- 
servation of NGC 1068 would be sufficient to constrain the 
polarisation angle with a statistical uncertainty of 5° if the 
polarisation degree is 10 per cent across 6 keV to 35 keV. For 
higher values of P, the precision on the measurement rises; a 
polarisation of 20 per cent constrains the angle with an un- 
certainty of 2.5°. At lower photon energies, the performance 
of the X-ray polarimeter is expected to be even better. The 
soft X-ray polarisation obtained by our modelling is at least 
comparable to and often significantly higher than 20 per 
cent as shown in Tabled Therefore, a future mid-size X-ray 
observatory with broad band polarimetry capabilities could 
directly measure the orientation of the innermost ionisation 
cones as described in Sect. 15. Il 

If we assume a vi ewing angle of i = 70° towards 
NGC 1068, as found by iHonig et all l|2007l ') from applying 
a clumpy torus radiative transfer model to infrared data, 
then the expected values for At/i given in Table [2] are also 



within the measurable limits. Should the inclination angle 
be hi gher as might be suggested by water maser observa- 
tions ^GammoreleEaDIliMlGreen^ then the 
possibility to measure A^ depends on the model details. We 
consider the values given in Table [2] as lower limits that are 
based on a conservative modelling approach. 

Note that if one decides to include both, a soft and 
a hard X-ray polarimeter (up to about 25 keV) in the 
large mission design of the International X-ray Observatory 
(jSarcons et al.ll201ll ). then measurements of the broad band 
polarisation still further improve. The method described 
here might then be applied systematically to a larger sample 
of objects. In summary, we therefore think that future X-ray 
polarimetry over a broad energy range is a promising tool 
to disentangle the details of the inner AGN geometry. 
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